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Abstract 

We consider an isotropic two dimensional harmonic oscillator with arbitrarily time- 
dependent mass M{t) and frequency Q(t) in an arbitrarily time-dependent magnetic 
field B{t). We determine two commuting invariant observables (in the sense of Lewis 
and Riesenfeld) L, I in terms of some solution of an auxiliary ordinary differential equa- 
tion and an orthonormal basis of the Hilbert space consisting of joint eigenvectors ipx 
of L, I. We then determine time-dependent phases a\{t) such that the 'ip\{t) = e^°'^ip\ 
are solutions of the time-dependent Schrodinger equation and make up an orthonor- 
mal basis of the Hilbert space. These results apply, in particular to a two dimensional 
Landau problem with time-dependent M, B, which is obtained from the above just by 
setting Q{t) = 0. By a mere redefinition of the parameters, these results can be applied 
also to the analogous models on the canonical non-commutative plane. 
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1 Introduction 



The study of time-dependent quantum problems has been of great interest in the hter- 
ature since the work of Lewis and Riesenfeld [HIS], in the late 1960s. They developed 
a completely general approach to the evolution generated by a time-dependent Hamil- 
tonian H{t); one re-obtains the results of the adiabatic approximation in the case of 
slowly varying H{t), and the results of the "sudden" approximation in the case of fast 
or even discontinuously varying H{t). In Refs.[l] and [2], the method was applied to 
a one-dimensional harmonic oscillator with a t-dependent frequency and to a charged 
particle in a t-dependent, axially symmetric electromagnetic field in three dimensions. 
For zero electric field, the latter reduces to a Landau problem with time-dependent 
magnetic field. Several other applications of the method have been worked out later; 
we just mention the application in [3\ to the anisotropic three-dimensional harmonic 
oscillator, the ones to the damped harmonic oscillator and to a [5] linear (in space) 
potential in one dimension. In the present work, we apply the approach of Lewis 
and Riesenfeld to a two-dimensional isotropic harmonic oscillator with time-dependent 
mass M and frequency Qina time-dependent (but space-independent magnetic field) B 
[formula ([3])] [the term Q'^x-x in the Hamiltonian may contain, in particular, an electric 
potential of the form (p{x) = ^jj^^^-x [B], [TB]] for = 0, this reduces to a Landau 
problem with time-dependent mass M, B. In the Coulomb gauge, the Hamiltonian 
can be rewritten as the sum of the Hamiltonian of a harmonic oscillator with suitable 
time-dependent mass ii and frequency oj and of the angular momentum L multiplied 
by a suitable time-dependent coefficient v [formula ([5])]. In either case the model can 
be immediately embedded in a three-dimensional one where B is directed along the 
third direction, and the motion in the direction of the latter is free; the latter would 
describe, in particular, a charged particle in a suitable t-dependent axially symmetric 
electromagnetic field (as in [2J, but with a possibly t-dependent mass). 

The results can be applied also to a two dimensional isotropic harmonic oscillator 
with time-dependent M,^},B on the (canonical) non- commutative plane [formula ([5])]. 
In fact, the Hamiltonian can be again rewritten in the same form ^ with new parame- 
ters w, v involving also the deformation parameter 9 of the non-commutative plane in 
their definition. Non-commutative space-time structure is an old subject dating back 
over fifty years to Snyder [8] . The analysis of field theories on non-commutative space- 
times has become a promising area in theoretical physics [9], and also non-commutative 
quantum mechanics has become subject for intense investigations. In the recent works 
[lOj the formulation, interpretation and applications of non-commutative quantum me- 
chanics have been investigated in configuration space, with the explicit example of the 
(time-independent) two-dimensional harmonic oscillator in [10], or general Hamiltoni- 
ans (including the one of the Landau problem) in fTT], where an approach to second 
quantization is also proposed. 
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In section 2, we define the model and use the socahed Bopp-shift to connect the 
noncommutative variables to the commutative ones. In section 3, we determine two 
commuting invariant observables L,I {L is the angular momentum) in terms of some 
solution of an auxiliary nonlinear ordinary differential equation. In section 4, we de- 
termine an orthonormal basis of the Hilbert space consisting of joint eigenvectors of 
L, /; the presence of the parameter of noncommutativity 6 does not modify much the 
results. Then, we determine suitable time-dependent phase factors such that their 
products with these joint eigenvectors make up an orthonormal basis of the Hilbert 
space consisting of solutions of the Schrodinger equation. Each phase factor can be 
split as the product of a dynamical one and geometric one (including the Berry phase) 
[?]. The latter includes, in general, the effects of the socalled adiabatic operator; the 
latter multiplied by the socalled dynamical operator and the non-adiabatic one gives 
a decomposition of the evolution operator, as explained e.g., in the introduction of [7] 
(see also [TS]). 

In the rest of this introduction, we briefly recall the general Lewis-Riesenfeld ap- 
proach to time-dependent quantum mechanics [2j. An operator I{t) is said to be an 
invariant if 

/(t)^|/(t) + l[/(t),/f(t)]=0 (1) 

(note that H itself is not an invariant unless it is time-independent). Here and in the 
sequel we use the Schrodinger [not the Heisenberg) picture. If I{t) is an invariant and 
'(/'(i) solves the Schrodinger equation 

in^-±{t) = H{t)m, (2) 

then also I{t)ip{t) does. If in addition I{t) is hermitian, then there exists a (in general 
time-dependent) orthonormal basis {ipx^n} of eigenvectors of / with real eigenvalues 
A S A; these eigenvectors are parameterized by A G A and if necessary by some 
additional label k £ K. Deriving the relation /(/Ja,k = ^fx,K with respect to t and using 
([1]), ([2]) one finds that all eigenvalues A are time- independent. Moreover, adjusting the 
phase, the ipx,K can be chosen so as to be a complete set of solutions of ([2]), provided 
/ does not contain the operator dt (for instance, operators of the form / = adt are 
not allowed). If H is time-independent and we choose a time-independent /, then eq. 
([1]) amounts to [/, H] = 0, and the above scheme amounts to finding a complete set of 
solutions of ^ in the form of eigenvectors of both I and H, as usual. 
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2 The model 



The Hamiltonian for the two dimensional isotropic harmonic osciUator in a magnetic 
field reads 



Hhm{t) 



2M{t) 
1 

2M{t) 



\p-eAit)] ■ [p-eA{t)] + ^^^^^^ x ■ x 



(3) 



B{t)x2Y + ip2-^Bit)xi 



we assume piecewise continuously time-dependent mass M, frequency and but mag- 
netic field B (assumed space-independent). The second line holds after expressing 
the vector potential in Coulomb gauge. As said, for O = this reduces to a Landau 
problem with time-dependent mass M, B, 

{p-eA{t)] • [p-eA{t)] 1 



HL{t) = 



2M{t) 



2M{t) 



Pi + ^B{t)x2Y + (p2-^B{t)xi^'^ 



■ (4) 



Wc use units such that h = 1 and Xa,Pa,^ arc dimcnsionlcss (this can be achieved 
renormalizing as usual the latter variables by dimcnsionful factors built out of some 
constant mass mo and frequency wq, and h). Using the canonical commutation relations 
[xajXb] = = [poiPb]) [a^ojPfe] = iSab it is straightforwad to check that Hhm{i) reduces 
to 

Hit) = 



1 Lo'it) iy{t) ^ 
p ■ p -\ — X ■ X —L 



1 



2 

2 , _2^ 



{pi + ^ {xi + xi) - u{t)L, 



where L = xip2 — X2P1 is the angular momentum, and 

li{t) := M(i), u'^{t) := M(i)J)2(t).i [^'^(*)]^ 



u{t) -r- 



eB{t) 



(5) 



(6) 



AM{t) ' 2M{t) ■ 

Our two-dimensional non-commutative space is defined by noncommuting coordi- 
nates Xi , X2 and associated momentum components P\ , P2 satisfying the following 
commutation relations: 



Xl,X2 


= ie, 




= h 


X2, P2 


= h 


Pl,P2 



0. 



(7) 



All these operators are hermitian. On such a non-commutative space we consider as 
a Hamiltonian that of a two dimensional isotropic harmonic oscillator in a magnetic 
field with f-dependent M,rt,B, 



+ 



M{t)n^{t) 



(8) 
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The aim is to find two independent invariants and then solve the Schrodinger equation 
within each eigenspace of the latter. We redefine the operators as follows 



Xi=Xi + -P2, X2 = X2--Pl, Pl = Pu P2=P2; (9) 

the associated commutations relations are 

[xi,X2]=0, [xi,pi]=i, [x2,P2]='i, \pi,P2]=0. (10) 

Such a redefinition is sometimes also called the Bopp shift. In terms of these new 
operators Hhmi't) reduces again to ([5]), where now 

^ 16M 4 ' 4M ' 8M 2 

(11) 

(we have omitted the time dependence for brevity). Note that in general ^^uj are 
t- independent only if all of M, B, are. 

Therefore for B{t) = 0, we obtain the harmonic oscillator with time-dependent 
frequency and mass, while for il(t) = we obtain the Landau problem with time- 
dependent frequency and mass, both on the commutative and noncommutative plane. 

An operator invariant of the model ([5]) is an operator that satisfies eq. ([TJ. 

3 Derivation of the operator invariants 

As x-x, p-p, x-p, p-x, and therefore H, are explicitly symmetric under rotations, they 
commute with the generator of rotations L, 

[L,x-x]=0, [L,p-p] = 0, [L,x-p]=0, [L,p-x]=0, [L,H]=0. (12) 

Since L is also time-independent, it fulfills ([1]) and is a first operator invariant. In order 
to find another one I{t) commuting with L, we make the Ansatz 

I(t) = a{t)x-x + l5{t)p-p + ^{t){x-p + p-x) + 5{t)L 

= a{t){xl + xl)+ p{t){pl+pl)+^{t){{xi,pi] + {x2,P2]) + mL, (13) 

where {xi,pi} = xipi+pixi, {x2-,P2} = X2P2+P2X2- As we need I{t) to be hermitean, 
= /, the coefficients a, /3, 7, 5 should be real. As said, 

[/,L] = 0, (14) 
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because / is manifestly symmetric under rotations. Using (jl2p and 

[xI,pI] =2i{xi,pi}, [xj,{xi,pi}] = Aixj, = -4ipf, 

[xI,pI] ='2i{x2,P2}, [xl,{x2,P2}] = ^ixl, [pI,{x2,P2}] = -^ipl, 

Eq. ([1]) reduces to 



d-2a;S = 0, (15) 

^J' 

y — I 



/3+-7 = 0, (16) 
7 + ia-a;2/3 = 0, (17) 



6 = 0. (18) 

Eq. (jlSp implies that (5 is a constant. The solution of the system (jl5p - (|18p with 5 = 1 
and all the other coefficients vanishing is the already determined invariant L. In the 
rest of the text, we choose 6 = for the other invariant. From (jl5p - (jl7p it follows 

l{j^-a(3)=0, (19) 

whence 

^2 _ Q,/3 = (20) 

where k is either a real or an imaginary constant (so that is real). 

Since I{t) and L commute, they have joint eigenvectors. One can find all of the 
latter using the same operator technique as in Ref. [1], that is the Dirac method of 
diagonalizing the Hamiltonian of a constant-frequency harmonic oscillator. To this 
end, we choose a positive k, stick to positive-definite solutions /3(t) and introduce an 
auxiliary real-valued function a{t) such that 

m = cy\t); (21) 

then by Eqs. ([16]) and ^ 

7 = —iJ,aa, a = & fj, H — (22) 

Replacing ([2T|) . (p2|) in (fTT]) . we find that the system ()15m7p is equivalent to the 
auxiliary equation 

a + ^a + —a = ^^. 23 
In terms of a real solution a{t) of (j23p the operator invariant (jlSp takes the form 

/(t) = ((T/iXl - Crpi) H ^X? + (cT/iX2 - crp2) H — ^T^^i. (24) 
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4 Eigenstates of I{t) and the phases 

We can introduce "raising" and "lowering " operators ai, 02, a\, 0,2 by 

aa{t) = ^= (c^l^Xa - apa + i-Xa) , = ^= (cT/J^Xa - OPa - i-Xc^ , (25) 

v2k ^ a / ■s/2k ^ a / 

with a = 1, 2. They fulfill the commutation relations 

[ai, a|] = 1, [02, 4] = 1, [ai, 02] = 0, [ai, 4] = 0. (26) 

Then the operator invariants /, L can be written as 

/_(t) = 2k ^4^*1 + 4^2 + 1^ , L = z(4ai — 4'^2)- (27) 

Instead of working with the , aJi it is convenient to work with the left and right 
circular annihilation operators A_ , Aj^ , that are respectively defined by 

A_(i) = -^(ai -^2), A+(i) = -^(ai+m2), (28) 

and the left and right creation operators, that are the hermitean conjugates of the 
above: 

Al (t) = i= (4 + ial), ^+ = 71 H - ^"2) ■ (29) 

The only non-zero commutators among A^, A^, A^_^, A[_ are then given by 

[A.,A^_] = 1, [yl+,4] = l. (30) 

Using the inverse relations 

ai = ^{A. + A^), a\ = %Al + Al), 



02 



2 V ■ 12 



we find 



and 



[I.,A±] = -2k, [I-,AI]=2k, 
[L,A±]=±A±, [L,Ai]=^Ai, 



(31) 

I-{t) = 2k{aIa+ + aIa_ + 1), L= {aIa_ - aIa+). (32) 



Let us assume ||0, 0) is a normalized state annihilated by ai, 02, or equivalently by A±: 

A_||0,0)=0, A+||0,0)=0. (33) 

For any r+, r_ G Nq (Nq stands for the set of nonnegative integers) we find 

I^{t)iAlY+{AlY- ||0, 0) = 2K(r+ + r_ + 1){aIY+{aIY- ||0, 0), 
L{t){AlY+{AlY-\\0,0) = (r+ -r_)(Al)^'+(^^, )^'-||0,0). 



(34) 



The pairs (r+,r_) are in one-to-one correspondence with the pairs {n,m), with n,m 
defined by 

n = r-i- + r_ G Ng, m = r+ — r_ G {— n, 2 — n, n — 2, n}. (35) 

The inverse of ([35]) are 

n + m n — m , ^ 

r+ = G No, r = ^— G Nq. (36) 

An orthonormal basis of the Hilbert space of states of the system consists of 

1 1 



n, m 



) =AA(At_)5("+™)(At )|("-'^)||0,0), M= — == (37) 



ri+m I / n— 'm i 



which by (|34p are eigenvectors of I-{t) and L: 

m) = 2K(n + l)||n, m), -L||n, m) = m||n, m). (38) 

The action of ^±,^j- on the basis reads: 



/ n+m *t II \ / n + m 
Aj^\\n,m) = Y — ^ — ||n — l,m— 1), A^j^\\n,m) = W — ^ hi ||n+l, m+(I^!^) 



. 11 , n — m ,, ^ ^, .+ II , n — m ^ ^ 
A^\\n,m) = Y — - — — l,m+l), A_\\n,m) = y — hi \\n+l,m—{^) 

By the general theory |2], one can transform the basis {||n,m)} into an orthonor- 
mal basis {ipn,m} consisting of solutions of the equation ([2]) applying suitable phase 
transformations 

V'n.m = e'"™'"W||n,m) (41) 
By the orthonormality of the basis, the Schrodinger equation 

d 

i— (e*"™'"||n,m)) = iJe*""-™ ||n, m) 



reduces [2] to the following equations for the time-dependent coefficients 



Or 



dOrn,n , n I ■ ^ iir \ II \ 

= {n,m\\ \ i— - H ] \\n,m}. 



(42) 



One can thus split am,n into the sum am,n = ctfn,n+am n °f ^ geometric phase (including 
the Berry phase phenomenon for slowly varying parameters) [121 [13} I14j (see also (THj ) 
and a dynamical phase a'^ „ respectively fulfilling 



i{n, m\\ 



d_ 
dt 



\n, m] 



a. 



m,n 



-{n, m\\H\\n, m) 



(43) 



To compute the matrix elements {n,m\\H\\n, m), we express H in terms of A±, A^^. As 
a first step, 



Xi 



cr{ai — a\) 



(44) 



in terms of A\,A^,A_ we find 



-'tcrfi 



P^ = ^ [A--A^_,+A+-Ai ) - ^ [A^+A^^+A^+A[ 



X2 



an, 



a.-aI-a++aI 
a_-aI-A:+aI 



' 2(T 



A_+Ai-A.-Al 



(45) 



and 



xi+ix2 = ^(aI-A^ 



P^-^P2 = '^.[Al-A-)-^[A-+A 



Xl — tX2 



Pi+ip2 = ft (AI-A+ ) [A++A 



it 



(46) 



Replacing these formulae in ^ we find 



H = ^[a-'fi^ — 2 



+ u^a^ 1 (^L^~ + AiA+ + 1) - z^(^L A_ - 



it 



it 



it 



+ 



,2 2 2 



/2k \/2(t 



,2 2 2 

^ ^ • 2k 



/2k ^/2o- 



(47) 



aUI. 
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Only the operators in the first hne contribute to the matrix elements {n^m\\H\\n^m) 
the ones in the second line lead to vanishing terms. So we finally find 

^2 



(n, m||i?||n, m) 

Now we wish to evaluate 
d 



n+l 
2k 



I 2 2 
+ W cr 



mv. 



(n,m||^||n,m) =AA(n,m||^ ) ^ ("+™) (A^) 5 (""'") || 0,0) 

|0,0) is determined up to a time-independent phase. 

On one hand, deriving (j38p with n = m = with respect to t, we find 



(48) 



(49) 



9/ 

5t 



|0,0) = {2k -I) ( ^||0,0) ) , 



L(|||0,0) 



0. 



(50) 



As the left hand-side of ()50P i is different from zero, ^||0,0) cannot vanish, but must 
be annihilated by L. We determine it by realizing Oq, ||0,0) in configuration space: 



^2k 



a 



] Xa + iada 



exp 



K i(T^\ X-X 



(51) 

in fact Vo is annihilated by both ai,a2, therefore also by A^,A-^-, and is normalized 
w.r.t. the scalar product (-0, V'') = / d'^x ip{x)'tp' (x). Deriving (j51|) 9 we find 



dipt] 
dt 



(x) 



ia^i ifia i&'^n 2Ka\ x-x a 
a a c7^ 1 2 a 



By the Stokes theorem and some straightforward computations we obtain 

,2 



d X da{\il^o{x)\ X' 







a 



dx\'il:Q{x)\ x-x = — ldx\ipQ{x)\ = — 



a 



From the two previous equations we conclude 



(o,o|| 



d_ 

dt 



|0,0) 



dt 



d\ tpo ipo 



/ iafi ifia i&'^fi 2k& 



x-x a 



2k 



■■ , Aff^o" .2 
(jcr H a 



(52) 



On the other hand, 

^ (yi;^)|("+™)(^t_ )!("-'") 



+ 



n + m dA^_ . i(„+^_3) t a(n-m) 
2 dt ^ ^ 



(53) 



10 



but from ([29]) it follows 



dA 



1 da\ ^ i dal 



dt ^ dt ^2 9t ' 
and by equations ([25|) . 



~dt 



1 da\ i 



V2 dt V2 dt' 



OA 



dt 



dAl 



1 



2^ 

in 

1 
2^ 



afi + fia + i— 2 (xi + ix2) - &{pi + ip2) 



aa + ^-a']A^ 



.KO 



aa 2 2iKa 
aa H & -\ 



A, 



aid + ila + i-^ ) (xi - 1x2) - &{pi - ip2) 



aa + 



/icr 



d'\A\-\ad + ^-a' + 



2iKa 



(54) 



(55) 



Inserting (|52]) - (l55]) into (I49l) and using (|23)1 we finally find 



(n, mil — lln, m) = i iaaix + fiaa — a fi) 

at 9.^ ^ / 



2k 



2k VA*'^^ 



22 -2 

Lo a — a II 



(56) 



as the operators A-|_,j4_ give vanishing contributions to the scalar product. Replacing 
SD and ([56]) in (02]), the latter takes the form 



dt 



-(n + 1)- 



+ mv 



(57) 



(there is a partial cancellation of the geometric and dynamical phase derivatives), which 
is solved by 



an,m.{t) = -{n + 1) / dt' 



— + m I dt' i'{t'). 



(58) 



5 Conclusions 

We have derived a class of exact invariants for the time-dependent, isotropic two di- 
mensional harmonic oscillator in a magnetic field, both on commutative and noncom- 
mutative plane. This applies also to the time-dependent Landau problem (which is 
obtained choosing zero harmonic oscillator frequency J7(t)). The angular momentum 
L (absent in the one dimensional case treated in Refs. [Ij) is trivially another invari- 
ant. Using these two invariants, we have constructed and an orthonormal basis of 
the associated Hilbert space consisting of solutions of the time-dependent Schrodinger 
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equation. The presence of the noncommutative parameter or of a time-dependent 
mass does not change too much the results with respect to the commutative space with 
constant mass, but induces a first order derivative term in the auxiliary equation ()23p . 
that is a general form of the Ermakov-Pinney equation in Refs.[T7]. When we set 9 = 0, 
then the equation (j58p are the same as the results in [2], but in two dimensions. For 
some choices of ^(t), exact solutions of (|23p can be found and then the corresponding 
solutions of the Schrodinger equation in close form. It would be also good to calculate 
the transition amplitude connecting any initial state in the remote past to any final 
state in the remote future with constant u; we hope to report on this point elsewhere. 

Note added in proof. In the final stage of publication of the present paper Prof. 
Dodonov has brought to our attention the very interesting papers, fl91, [2Qj, where 
among other things exact invariants in the form of generalized creation and annihila- 
tion operators for a general N-dimensional harmonic oscillator were constructed. 
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